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ABSTRACT 

A  maximum  likelihood  estimator  is  derived  for  monopulse 
radar  signals  which  have  been  coherently  integrated.  Field  data 
from  the  Millstone  Hill  radar  is  reduced  using  the  estimator  and 
the  results  are  compared  with  known  values  derived  from  an 
accurately  computed  ephemeris. 
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I. 


INTRODUCTION 


Several  radars  now  exist  which  are  capable  of  tracking  ob¬ 
jects  in  deep  space.  With  the  completion  of  the  SPADATS  modifi¬ 
cations  to  ALTAIR  this  radar  will  also  have  that  capability.  Deep 
space  is  an  inprecisely  defined  term  but  it  can  be  described  rela¬ 
tive  to  a  given  radar.  For  any  given  radar  a  satellite  in  deep 
space  has  the  following  attributes: 

i)  Angle  rates  are  small  and  thus  the  period  is 
probably  12  hours  or  more. 

ii)  The  SNR  is  low  and  is  probably  less  than  0  dB  on 
a  single  pulse. 

Of  most  importance  is  the  low  SNR  which  renders  single  pulse 
detection  unreliable.  In  general,  only  coherent  integration  can  be 
used  to  improve  this  SNR  to  a  point  where  reasonable  detection  and 
false  alarm  probabilities  can  be  achieved.  Such  integration  has 
been  used  on  the  sum  signals  in  all  deep  space  tracking  radars  to 
provide  target  detection  and  range  estimates.  However,  reduction 
of  the  tracking  data  to  provide  accurate  estimates  of  the  orbital 
elements  of  a  tracked  satellite  also  requires  accurate  measure¬ 
ments  of  the  monopulse  angles.  Several  ad  hoc  averaging  methods 
have  been  used  on  individual  radars  with  varying  degrees  of 
success.  In  this  report  a  new  estimator  of  the  monopulse  angle 
is  derived.  It  is  based  on  a  method  of  maximum  likelihood  but  it 
has  its  origins  in  an  analog  system  known  as  a  Normalized  Angle 
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Error  Detector.  A  detection  algorithm  falls  out  naturally  from 
the  derivation.  Bounds  on  the  performance  of  the  estimator  are 
derived  in  the  form  of  a  Cramer-Rao  bound  and  this  provides  results 
which  are  to  be  intuitively  expected. 

As  a  test  of  the  algorithm,  various  tracks  of  deep  space  ob¬ 
jects  were  conducted  by  the  Millstone  radar.  Data  recorded  during 
these  tracks  was  reduced  using  the  derived  algorithms  for  esti¬ 
mating  the  monopulse  error.  Very  close  agreement  was  found  be¬ 
tween  the  computed  monopulse  error  and  the  errors  derived  from  the 
radar  pointing  and  an  accurate  ephemeris. 

II.  DESCRIPTION  OF  THE  PROBLEM 

Objects  at  extended  range  may  not  be  detected  with  a  single 
pulse  as  the  SNR  may  be  at  or  below  the  noise  level.  However,  pro¬ 
vided  that  the  ephemeris  is  reasonably  well  known  and  the  target 
is  stable,  several  pulses  may  be  coherently  integrated  to  enhance 
the  SNR  sufficiently  such  that  reliable  detection  is  possible.  The 
ephemeris  is  required  in  that  the  gathered  data  must  be  adjusted 
for  any  accelerations  over  the  span  of  the  data.  The  traditional 
method  of  integration  is  through  a  Discrete  Fourier  Transform  (DFT) 
which  provides  coherent  integrations  for  a  set  of  linear  phase 
slopes.  If  a  detectable  return  is  in  the  samples  then  the  phase 
slope  which  most  nearly  approximates  the  actual  phase  slope  will 
provide  a  maximum  integrated  signal.  Failure  to  properly  compen¬ 
sate  for  accelerations  will  result  in  a  residual  quadratic  phase 


between  samples.  This,  in  turn,  will  result  in  a  spreading  of  the 
spectrum  -  the  output  of  the  DFT  -  over  several  cells  or  phase 
slopes. 

Spreading  of  the  spectrum  may  also  be  caused  by  differential 
motion  of  several  radar  scatterers  on  a  satellite.  This  motion 
would  be  the  result  of  rotations  of  the  satellite.  In  the  extreme, 
this  motion  could  spread  the  spectral  energy  over  all  aoppler  cells 
but  in  the  cases  of  interest  we  would  expect  the  spread  to  be  small 
when  compared  to  the  total  spectral  extent. 

For  this  analysis  the  target  will  be  modelled  as  a  series  of 
radar  scatters  each  with  a  different  doppler.  Further,  the  initial 
phase  of  the  return  from  each  scatter  will  be  random  with  a  uniform 
distribution.  Over  the  extent  of  the  data  gathering,  the  RCS  of 
each  scatter  will  be  constant.  As  no  attempt  is  being  made  to 
derive  detection  statistics,  the  scatterers  could  be  considered  as 
having  a  constant  RCS. 

For  a  single  radar  pulse,  a  monopulse  radar  derives  three 
signals,  a  'sum'  and  an  elevation  and  azimuth  'difference'  signal. 
Over  some  restricted  range  of  monopulse  angles,  the  monopulse 
angle  -  either  azimuth  or  elevation  -  is  proportional  to  the  quo¬ 
tient  of  the  relevant  difference  signal  divided  by  the  sum  signal. 
This  condition  applies  when  the  radar  is  in  proper  alignment  such 
that,  in  the  noiseless  case,  the  difference  signals  are  either  in 
phase  or  180  degrees  out  of  phase  with  the  sum  signal.  However, 


in  the  case  of  returns  from  deep  space  satellites  both  the  sum 
and  difference  channels  will  be  masked  by  noise  and  the  single 
pulse  estimate  of  the  monopulse  angle  will  be  meaningless. 

Several  approaches  are  available  for  estimating  the  monopulse 
angles  but  in  this  note  a  method  similar  to  that  used  in  an  analog 
Normalized  Angle  Error  Detector  is  used.  Consider  the  signals 
shown  in  Fig.  1.  In  this  case  the  phase  of  the  noiseless  dif¬ 
ference  signal  has  been  increased  by  90  deg.  The  tangent  of  the 
angle  0  between  the  sum  and  the  sum  plus  j  multiplied  by  the 
difference  is  the  monopulse  ratio.  Thus  an  estimate  of  0  leads 
to  an  estimate  of  the  monopulse  ratio. 

III.  ANALYSIS 

A.  Signal  Model 

Two  signals  shown  in  Fig.  1  are  considered  as  received 
vectors : 

-1  =  tZo  +  jAo'  E1  +  ^A1',,,EN-1  +  ^an-1^ 

and: 

-2  =  [lo  "  jAo'  E1  "  jAl' *  * ’ ,EN-1  ”  jAN-l]  (2) 

where : 

N  is  the  total  number  of  pulses 

However,  these  signals  are  composed  of  returns  from  several 
scatters.  In  addition,  there  is  a  model  of  an  underlying  trans¬ 
mitted  signal.  For  this  analysis  a  constant  phase  difference 
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and  difference  signals 


between  pulses  for  each  scatterer  corresponding  to  a  constant 
doppler  is  chosen.  This  underlying  model  applies  to  both  z-^  and 
z2-  Consider  K  scatterers  with  the  complex  amplitude  of  the 
composite  signal  (sum  plus  or  minus  difference)  of  {aR}.  The 
vector  of  complex  amplitudes  is  then: 


This  vector  incorporates  the  amplitude  and  initial  phase. 

Each  scatterer  has  an  individual  doppler  which  result  n  a 
constant  phase  shift  from  pulse  to  pulse  of  {<!>„}.  The  ch«  ~er 

J\ 

istics  may  be  incorporated  into  the  matrix  S  where 


N  is  the  number  of  integrated  pulses 
n  is  in  the  pulse  index  { 0,1, . . . ,N-1} 
k  is  in  the  doppler  index  {0,1, . . .K-l} . 

The  noiseless  models  of  the  received  signals  are  then 

r1  =  aSej0  (5 

r_2  =  aSe  (6 

B.  Noise  Model 

Noise  in  the  sum  and  difference  channels  should  be  uncorre¬ 


lated.  There  is  a  component  of  sky  and  ground  noise  in  both 
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*  - 


r* 


« 


« 
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signals  but  it  is  minor  compared  with  the  amplifier  noise  which 


is  independently  generated  in  each  channel . 

samples  in  each  channel. 

Consider 

the  noise 

nl  =  nZ  +  ^nA 

(7) 

n2  "  nz  -  ^nA 

The  correlation  is  then 

(8) 

E  nxn*  »  E|nzl2  -  E|nJ2  +  jE 

Provided  that: 

Re  n_n. 

Z  A 

(9) 

E  nznA  =  0 

and: 

(10) 

E|nJ2  =  E|nJ2  =  a2 

the  correlation  will  be  zero. 

(ID 

C.  Maximization  of  the  Likelihood  Function 

The  likelihood  function  may  be  written  as 

=  - i"N/2  exp  ~  { 11  g-i  ~  S-i  II  2  +  11  g-2  ~  L2  H  2}  (12) 

( 2ira  ^ )  2  a 

where 

z_^  and  z_2  are  t^ie  vectors  of  composite  received  signals 
r-^  and  r^  are  the  signal  models  from  Eq .  (5)  &  (6)  . 
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A  maximum  likelihood  estimate  of  the  parameters  is  that 
which  maximizes  the  likelihood  function  which  in  the  case  of 
Eq.  (12)  is  equivalent  to  minimizing  the  exponent.  Thus,  the 
estimate  is 


L(a,i,d) 


Min 
a,ip  ,6 


Expanding  Eq.  (13)  gives: 


where : 


and: 


(13) 

(14) 

(15) 


*  denotes  the  conjugate  transpose. 


The  signal  model  may  be  modified  by  quantizing  the  doppler 
shifts  into  an  equally  spaced  set  with  intervals  of  ^  cycles.  In 
this  case  the  set  of  dopplers  is  completed  to  a  full  set  numbering 
N  but  some  of  these  may  correspond  to  zero  amplitude  signals. 
Similarly,  the  vector  of  complex  amplitudes  is  expanded  to  N  by 
the  possible  addition  of  zero  amplitudes.  The  elements  of  the 
matrix  S  are  then: 


n  _  1  j2Trnk/N 

S|c  ”  "ST  ‘ 

The  rows  (and  columns)  of  S  now  form  an  orthonormal  system  such 
that 


SS*  =  I 


(17) 
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where 


I  is  the  unit  matrix. 

Applying  this  property  to  Eq.  (14)  gives: 


||rm||2  =  a  SS*  a* 
-  Hall2 


*  _  Jb  ic 

z^r^  =  z^S*  a  e 


*  Q-je 


=  -*  e^6  *  (2( 

The  transform  of  the  received  signal  is  a  Discrete  Fourier 
Transform  (DFT)  of  that  signal  and  may  be  denoted  as: 


-1  *  -Is* 


—2  ~  —2^*  ’ 


Equation  (14)  then  becomes: 


L  -  LQ-2  Re  (ux  e"^9  +<u2  ej9)a*  +  2||a||2  .  (23) 

The  function  L  in  Eq.  (23)  may  be  minimized  if  the  argument 
(phase)  of  a  is  chosen  as  follows: 

arg  a  =  arg  e  ]9  +  uj2  e-10  .  (24) 


Substituting  the  result  in  Eq.  (23)  gives: 


L  -  L0-2|(£1  e  ■’0+w2  e-’6||a|T  +  2||a||2, 
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where : 


is  a  vector  of  the  amplitudes  of  the  complex 
vector  a. 


Noting  that: 


Hill2  =  |a||a|T 


(26) 


The  expression  in  Eq.  (25)  is  minimized  by  selecting: 


'« 


|a|  =  e"j0  +  ±2  e j6  |  (27) 

Thus 

L  -  Lq  -  jHuj  e-j9  +^2  ej9  ||2  (28) 

Expanding  the  second  term  of  Eq.  (28)  provides: 

+  ||w2||2  +u.1»*  e‘j26  (29) 

(and  thus  L  is  minimized)  if: 

(30) 

(31) 

Equation  (31)  provides  a  maximum  likelihood  estimate  of  the 
angle  8  from  which  can  be  computed  the  monopulse  ratio. 

D.  Doppler  Detection 

Not  all  Doppler  cells  within  the  spectrum  will  be  occupied 
by  targets.  Thus,  some  threshold  must  be  applied  against  the 


i«i e  je  +“2 ej0 H2  =  H “i N 2 


Finally,  this  term  is  maximized 


0  =  j  arg 


or 


tan  20  = 


Im 

Re  a)  2 
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amplitude  of  a  particular  Doppler  return.  An  estimate  of  the 
amplitude  is  available  (Eq.  (27))  but  this  is  dependent  upon  a 
knowledge  of  the  angle  6 . 

A  recursive  method  of  detection  is  available.  If  an  initial 
estimate  of  zero  is  made  for  0,  then  Eq.  (27)  reduces  to  the 
traditional  detection  of  the  integrated  sum  channel.  For  those 
cells  in  which  a  signal  is  present,  the  sum-difference  angle  6 
may  be  computed  using  Eq.  (31) .  If  a  further  refinement  is 
required,  this  value  of  6  may  be  used  in  Eq.  (27),  new  estimates 
of  the  amplitude  made  and  a  new  threshold  detection  made  leading 
to  a  new  set  of  Doppler  cells  and  a  new  estimate  of  6 . 

E.  Bounds  on  the  Performance  of  the  Estimators 

Estimates  of  several  parameters  are  made  in  the  process  of 
estimating  the  sum-difference  angle  6 .  Each  of  these  estimates 
has  a  fundamental  limit  on  precision  of  the  estimate.  Lower 
bounds  on  this  precision  of  the  joint  estimation  process  is  given 
by  Cramer-Rao  bound.  Applying  this  bound  provides, 

o2 

o2(«)  i  - 2-y  (32) 

Hill2 

where : 

2  _ 

<?n  =  noise  power 

||  a ||  =  signal  power  in  detected  Doppler  cells. 


This  expression  then  reduces  to: 


°2<9>  5  2T1NR  <33> 

which  is  a  common  expression  for  the  performance  limit  of  an  angle 
detector. 

IV.  EXPERIMENTAL  RESULTS 

The  Millstone  radar  operates  in  a  mode  where  the  theory 
derived  in  previous  sections  is  particularly  applicable.  Many 
of  the  satellites  tracked  by  that  radar  require  coherent  integra¬ 
tion  to  extract  a  signal  of  sufficient  SNR  for  tracking  and 
detection.  Data  from  that  radar  has  been  analyzed  and  in  the 
following  section  the  results  of  these  analyses  are  presented. 
These  results  appear  to  be  accurate  bearing  out  the  usefulness 
of  the  estimator. 

A.  Data  Gathering 

An  initial  track  was  made  of  an  object  to  establish  an 
accurate  ephemeris.  With  this  in  hand,  data  was  taken  in  three 
parts : 

i)  An  azimuth  scan  from  -0.3  degree  to  0.3  degree 
at  0.1  degree  steps  wit  i  zero  elevation  offset. 

ii)  An  elevation  scan  similar  to  the  azimuth  scan 
with  zero  azimuth  offset. 
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iii)  Two  corner  points  with  offsets  at  -0.2  degree 
of  azimuth  and  elevation  and  0.2  degree  of 


•a 

^  • 


azimuth  and  elevation. 

Data  in  excess  of  3500  points  was  taken  at  each  station. 

B.  Processing  Algorithms 

At  each  station  the  data  was  processed  with  the  algorithm 
derived  in  the  previous  section.  The  following  narrative 
describes  the  processing: 

i)  For  each  data  station  and  for  both  azimuth  and 
elevation  the  derived  samples  of  Sum  +  j  Difference 
and  Sum  -  j  Difference  were  formed. 

ii)  Partitions  of  512  and  1024  samples  of  both  signals 
were  formed  and  transformed  using  an  FFT  function. 

iii)  An  initial  estimate  of  a  zero  sum-difference  angle 
was  made.  The  amplitude  of  the  sum  of  the  two 
signals  was  formed  and  thresholded  at  10  and  15  dB 
above  a  previously  estimated  noise  level. 

iv)  For  those  spectral  cells  which  were  above  threshold 
an  estimate  of  the  sum-difference  angle  was  made 
using  Eq.  (31) . 

v)  Estimates  of  the  monopulse  angle  were  made  using 
the  Millstone  monopulse  coefficient  of: 

Monopulse  angle  =  0.4  tan©  (34) 


.3 
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c. 


Review  of  Experimental  Results 


The  processed  data  and  the  results  appeared  to  be  in  accord 
with  the  predictions.  Precision  of  the  estimates  is  clearly  de- 

» 

*  pendent  on  the  SNR  and  is  suitably  low.  A  biasing  of  the  estimates 

appears  to  be  present  but  this  too  is  relatively  small. 

1.  Azimuth  Scans:  Tables  1  and  2  show  the  estimates 
^  and  precision  for  the  azimuth  scan  detected  at  15  dB  above  the 

noise.  The  SNR  figure  represents  an  estimate  of  the  total  SNR  of 

the  whole  system,  that  is,  the  signal  power  divided  by  the  relevant 

m 

noise  power. 

It  should  be  noted  that  a  bias  exists  and  that  it  is  entirely 
reproduceable  with  the  two  transform  size. 

|  *  TABLE  1 

I  AZIMUTH  MONOPULSE  ANGLES:  1024  POINT  INTEGRATION 


True 

(deg) 

Estimate 

(deg) 

S.D. 

(deg) 

SNR 

(dB) 

-0.30 

-0.24 

0.02 

27.58 

-0.20 

-0.15 

0.01 

27.46 

-0.10 

-0.06 

0.01 

29.95 

0.00 

0.01 

0.01 

28.82 

0.10 

0.09 

0.  01 

29.84 

0.20 

0.19 

0.01 

28.90 

0.30 

0.27 

0.0 

23.81 

TABLE  2 

AZIMUTH  MONOPULSE  ANGLES:  512  POINT  INTEGRATION 


True 

(deg) 

Estimate 

(deg) 

S.D. 

(deg) 

SNR 

(dB) 

-0.30 

-0.24 

0.03 

25.05 

-0.20 

-0.15 

0.01 

24.53 

-0.10 

-0.06 

0.01 

24.34 

0.00 

0.01 

0.01 

26.67 

0.10 

0.09 

0.01 

23.79 

0.20 

0.19 

0.02 

23.23 

0.30 

0.28 

0.03 

23.37 

Figure  2  gives  a  plot  of  the  estimated  azimuth  monopulse 
angle.  The  bias  term  is  clearly  apparent.  However,  the  bias 
could  be  reduced  by  modifying  the  slope  which  is  equivalent  to 
recalibrating  the  Millstone  monopulse  constant  of  0.4  deg/v/v/. 

2.  Elevation  Scans:  Results  very  similar  to  those 
for  the  azimuth  scan  were  obtained.  The  bias  term  is  small  and 
thus  the  recalibration  would  be  smaller.  It  would  also  be  in  the 
opposite  direction.  Tables  3  and  4  show  the  elevation  monopulse 
for  two  transform  sizes  and  a  15  dB  detection  level.  Figure  3 
is  a  plot  of  the  estimated  monopulse  angle. 
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Fig.  2.  Average  azimuth  monopulse  angles:  1024  point 
integration. 
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ESTIMATED  ANGLE  (deg.) 


-.40  -.20  0  .20 


TRUE  ANGLE  (deg.) 


Fig.  3.  Average  elevation  monopulse  angles:  1024  point 
integration. 
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TABLE  3 

ELEVATION  MONOPULSE  ANGLES:  1024  POINT  INTEGRATION 


True 

(deg) 

Estimate 

(deg) 

S.D. 

(deg) 

SNR 

(dB) 

-0.30 

-0.42 

0.00 

24.95 

-0.20 

-0.24 

0.01 

28.08 

-0.10 

-0.12 

0.00 

29.39 

0.00 

-0.01 

0.01 

28.82 

0.10 

0.10 

0.00 

29.12 

0.20 

0.22 

0.00 

25.51 

0.30 

0.3/ 

0.02 

25.24 

TABLE  4 

ELEVATION  MONOPULSE  ANGLES:  512  POINT  INTEGRATION 


True 

(deg) 

Estimate 

(deg) 

S.D. 

(deg) 

SNR 

(dB) 

-0.30 

-0.43 

0.03 

22.43 

-0.20 

-0.24 

0.02 

25.06 

-0.10 

-0.12 

0.00 

26.38 

0.00 

-0.01 

0.01 

26.67 

0.10 

0.10 

0.01 

26.48 

0.20 

0.21 

0.02 

24.92 

0.30 

0.38 

0.01 

21.86 

3.  Effects  of  Detection  Threshold:  Changes  in  the  de¬ 
tection  threshold  will  affect  the  precision  of  the  estimate.  How¬ 
ever,  the  total  effect  will  depend  upon  the  distribution  of  energy 
in  the  various  spectral  lines.  A  low  threshold  will  pick  up  low 
energy  spectral  lines  but  it  will  also  pick  up  noise.  Higher 
thresholds  will  restrict  the  number  of  noise  returns  but  may  also 
exclude  signal  energy. 
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A  sample  spectrum  of  a  sum  signal  is  shown  in  Figure  4. 

Energy  is  largely  restricted  to  a  small  number  of  cells  with  a 
peak  which  is  approximately  25  dB  above  the  average  noise  floor. 

For  the  case  examined  (10  dB  and  15  dB  thresholds)  very  little 
difference  may  be  noticed.  The  same  biases  exist  and  are  repro- 
duceable.  Table  5  shows  the  estimates  for  both  azimuth  and  ele¬ 
vation  channels. 


TABLE  5 

TRUE  AND  ESTIMATED  MONOPULSE  ANGLES  FOR  1024  POINT  INTEGRATION 
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SPECTRAL  POWER  (dB) 


V. 


SUMMARY 


An  estimator  has  been  devised  which  combines  multiple 
coherent  radar  returns  (both  sum  and  difference)  to  produce  an 
estimator  of  the  monopulse  angle.  The  estimator  is  derived 
through  maximum  likelihood  principles  rendering  it  optimum  in 
that  sense.  It  is  easily  implemented  using  standard  processing 
functions  of  FFT's  and  complex  arithmetic. 

Experimental  data  taken  by  the  Millstone  radar  and  when 
processed  according  to  the  derived  algorithms  has  provided  results 
which  are  reasonably  accurate.  However,  they  point  to  biases 
which  are  probably  within  the  radar  system  but  which  could  be 
removed  by  additional  calibration. 
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